Polymer chain stiffness versus excluded volume: A Monte Carlo study of the 
crossover towards the wormlike chain model 
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When the local intrinsic stiffness of a polymer chain varies over a wide range, one can observe both 
a crossover from rigid-rod-like behavior to (almost) Gaussian random coils and a further crossover to- 
wards self-avoiding walks in good solvents. Using the pruned-enriched Rosenbluth method (PERM) 
to study self-avoiding walks of up to Nb = 50000 steps and variable flexibility, the applicability of 
the Kratky-Porod model is tested. Evidence for non-exponential decay of the bond-orientational 
correlations (cos 9(s)) for large distances s along the chain contour is presented, irrespective of chain 
stiffness. For bottle-brush polymers on the other hand, where experimentally stiffness is varied via 
the length of side-chains, it is shown that these cylindrical brushes (with flexible backbones) are not 
described by the Kratky-Porod wormlike chain model, since their persistence length is (roughly) 
proportional to their cross-sectional radius, for all conditions of practical interest. 
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I. INTRODUCTION AND MOTIVATION 



One of the most basic characteristics of macro- 
molecules with linear "chemical architecture" is chain 
flexibility (or lack thereof, stiffness) [H-Q- While many 
synthetic polymers are fully flexible, and the statistical 
properties of their conformations under good solvent con- 
ditions have been extensively investigated [D-Q , recently 
also semiflexible polymers have found much interest, in 
particular since important biopolymers such as DNA, 
some proteins, rodlike viruses, or actin filaments belong 
to this class Moreover, also synthetic polymers 

such as polyethylene exhibit some stiffness over short dis- 
tances along the chain. Particular interest in the problem 
of chain stiffness has arisen due to the discovery that in 
macromolecules with "bottle-brush architecture" 
their stiffness can be varied over a wide range by changing 
the length and grafting density of side chains. Since short 
bottle-brush polymers in solution may exhibit liquid crys- 
talline type ordering [l3j , and their structure is very sen- 
sitive to various external stimuli, various applications for 
such molecules have been proposed (from building blocks 
of nanostructures to actuators and sensors etc. [Til [To]]). 
Also in a biological context semiflexible macromolecules 
with bottle-brush architecture are found [16], and have 
interesting functions such as lubrication in mammalian 
joints 

While for fully flexible polymers under good solvent 
conditions excluded volume interactions (loosely speak- 
ing, monomers of a polymer chain cannot "sit" on top of 
each other) are of central importance [lJ-Q , the standard 
model to describe semiflexible chains ignores them com- 
pletely (18l422i |. This "standard model" is the "wormlike 
chain model" of Kratky and Porod [l8|, described by a 
Hamiltonian (in the continuum limit) 



H 



dt 



(d 2 r(t)\ 
V dt 2 J 



(1) 



where the curve r(t) describes the linear polymer of con- 
tour length L, and the parameter k describing the bend- 
ing stiffness is related to the "persistence" length £ p as 
k = tphsT '. For chain molecules in the absence of ex- 
cluded volume, where distances between monomers that 
are far apart satisfy Gaussian statistics, one can show 
that the orientational correlation function between bond 
vectors decays exponentially with s = t/£ b 



(Si ■ a i+s ) = ilexp(-s£ b /£ p ) 



(2) 



where £ b = |oi| is bond length between two subsequent 
monomers along the chain, <Zj = fj — fj_i. Note that s 
is dimensionless and for a model with discrete monomers 
just denotes the difference in the labels i, i + s of the 
monomers. The Kratky-Porod model then yields for the 
end-to-end vector R e of the semiflexible chain [l8| (L = 
N b lb) 



(R 2 e )=2£ p L{l- £ j-[l-eM-L/£ P )}} 

which shows the standard Gaussian behavior ((R 2 ) = 
2£ p L = 2£ p £ b N b ) for L -> oo, while for L « £ p the 
rod- like behavior (R 2 ) = L 2 results. 

However, for good solvent conditions and long enough 
chains cqs. ©, and (J3|) cannot be correct: rather we must 
have [H, for N b -)• oo, 

(oj • a i+s ) oc s'P , = 2 - 2v w 0.824 , s* < s < N b 

(4) 



(3) 
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and rim m m 



(Ri) = 2t p , R t b Ni 



2// 



(5) 



(a) 



where v w 0.588 [H| is the exponent describing the 
swelling of the chains due to excluded volume. If s is 
no longer much smaller than N b , the power law eq. (j4|), 
gradually crosses over to a faster decay that depends on 
Nb In eq. ([5]), we have written a tentative gen- 

eralization of the Kratky-Porod result to introduce an- 
other, effective persistence length £„ r [25[. The question 
that we ask in this paper hence is, how can one recon- 
cile eq. ([2]) with eq. as well as eq. (|3j) with eq. ([5])? 
Tentatively, one might expect that for semifiexible chains 
eq. @ still holds up to some characteristic, large value 
s*, where then the crossover to eq. (g]) takes place [26j j ; 
but such a hypothesis is by no means evident, and re- 
mains to be proven, and the value of s* remains to be 
predicted. Similarly, one might expect that eq. ([3]) holds 
only for Nb up to some value N£, and then a crossover 
to eq. ([5]) takes place. Using a Flory [27[ argument, Netz 
and Andelman [28[ suggested that this is the case and 
(implying v = 3/5) 



(b) 



0.1 



0.01 





=0.1 


\\ Pb= 


=0.2 


\\ Pb= 


=0.4 




= 1.0 






V **.* \ x 




V \ A \ 

*• V V * A X 

\ 




X Xv 

V w ~ **l(v,/„ 



0.1 



0.01 



0.001 



10 



15 

Slh 



20 25 



30 




10 



100 



(6) 



As expected, this result for I v ,r is equivalent to the clas- 
sical result (H-ill {R 2 e ) 1/2 oc 4(^ P /4) 1/5 iV b 3/5 for semi- 
flexible macromolecules in the limit Nb — > oo. 

In the present work, we hence wish to check whether 
such ideas about these crossovers apply, and if so, study 
their detailed behavior. Finally, we shall discuss whether 
or not these ideas have some bearing on the problem of 
the persistence length of bottle-brush polymers 25J . 



II. MODEL AND SIMULATION RESULTS 

We carried out Monte Carlo simulations of self- 
avoiding walks (SAWs) on the simple cubic lattice, ap- 
plying an energy £(,(1 — cos 9) if a bond orientation differs 
by an angle relative to the preceding bond (of course, 
on our lattice only 8 = or 9 = ±ir/2 is possible). 
Using the pruned-enriched Rosenbluth method (PERM) 
[32T , the partition function of the SAWs of Nb steps 
with -/Vbend local bends (where 9 — ±7r/2) is written as 
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Z Nb ,N hen M) = E C(N b ,N hend )q^ (7) 

config 

where qb = exp[— Sb/ksT] is the appropriate Boltzmann 
factor (qb = 1 for ordinary SAWs). We obtained data 
for C(N b , N hcnd ) for N b up to N b = 50000, and vary q b 
over a wide range as well, 0.005< q b < 1.0. In addition 
we have continued our simulations (25j j of bottle-brush 
polymers using the bond fluctuation model. 



FIG. 1. Semi-log plot of the bond vector correlation function 
(cos 8(s)) vs. the chemical distance IbS along the chain (here 
it is the lattice spacing which is our unit of length, It = 1) 
for 0.1 < qt < 1.0 (a), and for 0.005 < qt < 0.05 (c). Part (b) 
shows the same data as part (a) but on a log-log plot. The 
straight line indicates a fit of eq. ^ to the data for qt — 1 
while straight lines in (a) (c) are fits of the initial exponential 
decay to eq. ([2]). Part (d) shows the distribution function 
P(n a tr) of straight pieces of the chain without kinks, together 
with fits of P(ristr) oc exp(— n str /n p ). Nb = 50000 was used 
throughout. 
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Fig. Q] shows that two distinct patterns of behav- 
ior emerge: for flexible or only moderately stiff chains, 
0.1 < qb < 1-0, eq. ([2]) has hardly any significant regime 
of applicability, there are just a few discrete values for 
small s(s — 1, 2, 3, • • • ) that can be fitted by eq. ([2]), and 
then pronounced deviations from a simple exponential 
behavior occur (fig. QJi). On the other hand, for large 
s (s > 20) the power law {eq. (g])} provides a good fit 
(fig- lib)- But it is also evident that the deviations from 
the power law for small s become the more pronounced 
the smaller qb is. Of course, the crossover value s* is not 
sharply defined, but rather the crossover is smeared out 
over some range in s. For small qb (0.005 < qb < 0.05) 
the initial exponential decay {eq. ((2|)} becomes better 
visible, fig. [TJ;, and chain lengths even larger than stud- 
ied by us would be needed to still clearly identify the 
power law for large enough s. While the exponential de- 
cay of the bond correlations {eq. ((2])} works only for a re- 
stricted range of s(s <C s*) and only for very stiff chains 
(qb <C 1), a simple exponential decay always is found 
for the probability P(n str ) that sequences of n str subse- 
quent bonds without kink occur, for large n str (fig. [TJl, 
P(n s tr) oc exp(— n stI /n p ) where n p then is defined as de- 
cay constant. For qb <C 1 we find that £ p = n p £b, within 
numerical error, and furthermore £ p oc q^ 1 holds then. 

The chain linear dimensions tell a similar story (fig. [2]) : 
for flexible chains, qb = 1.0 and qb = 0.4, there is a mono- 
tonic increase towards the plateau that yields the value 
£p,R ( ec l- ©), fig- [2^. However, the resulting values for 
£ Pi b, are small, unlike the case of bottle-brushes (fig.EJ)), 
with one side chain per backbone monomer, where with 
increasing side chain length N rather large values of £ Pi r 
can be achieved 25] . The strong initial rise seen for small 
qb in fig. [2^,, however, is indeed compelling evidence that 
for small qb a rod-like behavior is found. But also for large 
qb in fig. [2ji, and also in the case of the bottle-brushes 
(fig. EJd), one observes for small Nb a strong increase of 
(Rg) with Nb, but not as strong as it would occur for 
rods. 

Fig. |3] now tests the applicability of the Kratky-Porod 
formula, eq. (O, to these data; extracting £ p from the 
exponential fit of eq. ([2]) to the data in fig. QJ, there is 
no adjustable parameter whatsoever! Here the data for 
qb > 0.2 were omitted, since they clearly do not show 
any trace of the Gaussian behavior implied by eq. ([3]). 
For small Nb and small qb the success of eq. ([3]) indeed 
is remarkable (fig.EK), however, for large Nb we see that 
the data do not really settle down at the "plateau" value 
implied by the Gaussian behavior, but rather start to 
rise again, the crossover to the excluded volume behavior 
sets in. Again, this crossover is spread out over about 
a decade in Nb (and for qb = 0.005 it starts only for 
N b > 50000!). 

In experiments, information on chain stiffness of poly- 
mers, in principle, is accessible from an analysis of the 
static structure factor of single chains where q 

is the scattering vector. From a "Kratky plot", qS(q) 
vs. q, one finds a position where S(q) is maximal, 
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FIG. 2. Rescaled mean square end-to-end distance 
{R^} /(2£bN^) plotted against the chain length Nb for semi- 
flexible chains with £ b — 1 and variable qb (and hence £ p ) 
on a log- log plot (a), and analogous quantity for the back- 
bone of a bottle-brush polymer [25|], using the bond fluctua- 
tion model with excluded volume interactions only which has 
lb = 2.7 lattice spacings, and various side chain lengths N, on 
a linear-linear plot (b). The slopes shown in (a) illustrate the 
exponents expected for rods ((Re) oc N b ) or Gaussian chains 
((Re) oc Nb), respectively. Lengths are measured in units of 
the lattice spacing throughout. 



9max oc l/i? gyr at small q, where -R gyr is the gyration ra- 
dius of the polymer, and a power law qS(q) oc 
for chains exhibiting excluded volume statistics at larger 
q. For very stiff chains, one would expect that this power 
law crosses over to the Gaussian behavior qS(q) oc q^ 1 at 
large q and ultimately a further crossover to the scatter- 
ing function from straight rods, qS(q) = const., should 
occur, fig. [4^, shows that the latter regime is only found 
for the stiffest chains that we could study (qb = 0.005), 
and the onset of this so-called "Holtzer plateau" [34[ in 
the Krat ky p lot is rather gradual: contrary to some sug- 
gestions [35[ HH, estimation of the value q* where this 
onset occurs is not suited as an accurate method to ex- 
tract l p (= 1/q* )■ Fig. @Ji also shows that the expected 
power laws cannot be clearly identified either, because 
the crossover between them is too gradual. Fig. [4Jd re- 
veals that also the position i7 max where qS(q) has its 
maximum does not (yet) exhibit the expected power 
laws [28| : the reason is that for small qb the chain length 
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FIG. 3. Plot of (Rl)/2l b N b vs. N b , on linear-linear scales, 
for Nt < 40 (a) and on log- log scales, up to N b = 50000 (b), 
using only the data for 0.005 < Oh < 0.1. Dotted curves refer 
to the discrete chain model of [2l| . 



Ni, = 50000 still is by far too short that the mean square 
radii (R^) (cf. fig. [5k) and (Rg yl ) can reach their asymp- 
totic behavior. The smaller q b becomes (and the larger 
hence £ p is) the more Gaussian character of the chains 
do we expect for fixed Nb = 50000, and hence we see 
a gradual crossover from <7 max oc £ p to gwx oc l p ^ 2 
with increasing £ p on the log-log plot. 

Returning now to the behavior of bottle-brush poly- 
mers (fig. [2}d), we emphasize that pronounced local stiff- 
ness occurs there by a completely different mechanism, 
namely as a consequence of the local thickening of the 
(coarse-grained) cylindrical chains. This thickening with 
increasing side chain length N is evident from the radial 
density distribution in the plane (locally [371 ] ) perpendic- 
ular to its backbone (fig. [5k). Normalizing p(r) in fig. [5k 

oo 

such that N = 2~k J p(r)rdr, the cross -sectional radius 
o 

Rcs(N) can be defined as i?^ s = 2tt J p(r)r 3 dr. We now 
view the bottle-brush polymer as a sequence of blobs with 
diameter 2R CS (N). Fig. \5jp shows the construction. On 
the ordinate the values for 2R CS (N) are indicated, and 
we can identify how many backbone monomers Shiob(N) 
occur per blob by requiring Ar(sbiob(A0) = 2R CS (N). 
We find Sbiob = 6, 10, 12 and 14, for N = 6, 12, 18 and 
24, respectively. If we rescale Nb with Sbiob, and (R^ b ) in 



FIG. 4. (a) Rescaled Kratky-plot, log-log plot of 
qS(q)/[qS(q)] m!1 .x vs. g/g m ax, for several choices of q b as in- 
dicated. Straight line indicates the effective exponent that 
occurs in the regime of intermediate wavenumbers q. The 
position q max where qS(q) exhibits its maximum [g5 , (q)] max , 
and the height of this maximum is plotted in (b) as function 
of £ p . Straight lines indicate the exponents g max oc £p 1//5 and 
Qmax oc l p that one expects according to Netz and Andel- 
man [2S| in the excluded volume regime and Gaussian regime, 
respectively. All data have been taken for iVi, = 50000. 



fig.rj) with 2£ b £ PiR N^, we find that all data of fig.rjj fall 
on a universal curve (fig. Hb)!. The scaling implies that 
for bottle-brush polymers the large values of the effective 
persistence lengths £ Pt n are entirely due to the side-chain 
induced thickening of these flexible cylindrical brushes, at 
least for the range of side chain lengths accessible in the 
simulations. But this range nicely coincides 38| with the 
range accessible in experiments 0, [ll], fl3l4l5l .l35l l3fj| . l39j . 

Our findings that the persistence length £ p of bottle- 
brush polymers is of the same order as the cross-sectional 
radius R cs agrees with the result of Birshtein et al. (ipj 
but disagrees with later scaling theories QilUlj]. In fact, 
if £ p oc R cs the second virial coefficient v-i scales as [28[ 
V2 oc £ p R cs oc £p, and the Flory theory would yield 

an cnd-to-cnd distance [H R e oc (u 2 /^) 1/5 ^ /5 ^ /5 = 



£p /5 £ 3 b /5 N b /5 . A Gaussian behavior, R\ = £ b £ P N b , is 
only expected to occur for £ p /£b < Nb < Nj*. How- 
ever, equating the above two expressions for R e yields 
N£ oc £ p /£b as well, i.e., in agreement with our findings, 
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FIG. 5. (a) Radial monomer density distribution p(r) in 
planes locally perpendicular to the backbone of bottle-brush 
polymers with backbone lengths Nt = 1027 plotted vs. dis- 
tance r for side chain lengths N = 24, 18, 12 and 6. (b) Spa- 
tial distance Ar(s) between monomers a chemical distance s 
apart plotted on a double-logarithmic scale. The construction 
Ar(s) = 2R CS (N) is shown, (c) Rescaled mean square end-to- 
end distance of bottle-brush backbones, {R\ b ) /(2£t£ p .RN^") 
plotted vs. the rescaled chain length iVb/sbiob as described in 
the text. Note that our length unit (lattice spacing) physically 
corresponds to about 0.3nm [38l ] 



preasymptotic Gaussian regime does not occur for bottle- 
brushes. The scaling theories on the contrary imply that 
R cs oc iV 3 / 4 and l v oc TV 15 / 8 , and hence clearly would 
predict the existence of a Gaussian regime. Furthermore, 
even for much larger side chains the behavior R cs oc 
could not be verified [Hj]. Although we clearly cannot 
exclude that the scaling theory might become valid for 
side chain lengths of the order of N « 10 3 or larger, 
this clearly is completely irrelevant for both experiments 

[1 [TJ EM! EE EE Hf and simulations J2I M, M, S3 ■ 

In fact, self-consistent field calculations [43j have shown 
that side chain interactions do have an important effect 
onto the apparent persistence length only if N > 1000. 



III. CONCLUSIONS 

In summary, we have verified two distinct scenarios 
for local stiffness of polymer chains: if chain rigidity is 
caused by an energy term that suppresses chain bending 
to a large extent, without increase of the cross-sectional 
radius of the chains, one finds a regime of wormlike chain 
behavior, described by the Kratky-Porod model {e.g., 
eq. ([3])} for chains up to a chain length N b (where N b 
can be estimated by a Flory argument, cf. eq. (HJ) 28]), 
where a gradual crossover to excluded- volume-dominated 
behavior {eq. ([5])} occurs. However, if chain rigidity is 
caused by an increase of local chain thickness (as oc- 
curs for bottle-brush polymers) , no preasymptotic Gaus- 
sian regime compatible with the Kratky-Porod model oc- 
curs. However, another preasymptotic regime occurs, for 
2 < AV Sblob < 20, where (R 2 e 
from the rod-like behavior (B% 
volume behavior, (R^ b ) oc (fig. [5b)- In a regime 

where the contour length exceeds the persistence length 
only a few times, the Kratky-Porod model can be fit- 
ted to experimental data, excluded volume not yet being 
very important. However, other experiments may fall 
in a reg ime where excluded volume effects are impor- 
tant [44H46| . In [H, [iH the importance of excluded vol- 
ume was inferred from the concentration dependence in 
semi-diilute solution. Clearly, it is important to carefully 
examine to which regime experimental data belong. 
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